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ABSTRACT 

To date, the only far-infrared spectroscopic observations of ultralu- 
minous infrared galaxies have been obtained with the European Space 
Agency’s Infrared Space Observatory Long Wavelength Spectrometer. 
The spectra of these galaxies are characterized by molecular absorption 
lines and weak emission lines from photodissociation regions (PDRs), but 
no far-infrared (A > 40 pm) lines from ionized regions have been detected. 
ESA’s Herschel Space Observatory, slated for launch in 2007, will likely 
be able to detect these lines in samples of local and moderate redshift 
ultraluminous galaxies and to enable measurement of the ionization pa- 
rameters, the slope of the ionizing continuum, and densities present in the 
ionized regions of these galaxies. The higher spatial resolution of proposed 
observatories discussed in this workshop will enable isolation of the cen- 
tral regions of local galaxies and detection of these lines in high-redshift 
galaxies for study of the evolution of galaxies. Here we discuss evidence 
for the effects of absorption by dust within ionized regions and present the 
spectroscopic signatures predicted by photoionization modeling of dust- 
bounded regions. 
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1. Background 

Prior to the launch of the European Space Agency’s Infrared Space Observatory 
(ISO) in 1995, Voit (1992) showed how mid- and far-infrared fine-structure lines could 
be used to constrain the electron densities, extinction, and shape and ionization pa- 
rameters of the central ionizing sources in ultraluminous infrared galaxies (ULIRGs). 
Moreover, the ground-based work of Roche et al. (1991) showed that the mid-IR spec- 
tra of the nuclei of galaxies could be placed into three classes: those with aromatic 
feature emission, featureless, and those with silicate absorption, typically associated 
with optically identified H II region, Seyfert 1, and Seyfert 2 nuclei respectively. 
Building on this early work on optically selected starburst and AGN galaxies, Genzel 
et al. (1998) constructed a diagnostic diagram of the ratio of high-to-low ionization 
fine-structure lines vs. the equivalent width of the 7.7 gm aromatic feature emission 
based on which they concluded that 70 - 80% of ULIRGs are powered predominantly 
by starbursts and 20 - 30% are powered by a central AGN. They attributed the weak- 
ness of the mid-infrared fine-structure lines relative to the infrared luminosity to the 
effects of extinction. 

Far-infrared spectroscopy of a small sample of IR-bright and ultraluminous galax- 
ies taken with the ISO Long Wavelength Spectrometer (LWS) has revealed a dramatic 
progression extending from strong fine-structure line emission from photoionized and 
photodissociated gas in the starburst galaxy Arp 299 (Satyapal et al. 2002) to faint 
[C II] 158 pm line emission and absorption in lines of OH, H 2 0, CH, and [O I] in 
the ULIRG Arp 220 (Figure 1; Fischer et al. 1999). With the progression towards 
weak emission line strengths, no trend in density or far-infrared differential extinc- 
tion is indicated (Figure 2), i.e. the temperature-insensitive [O III] 52 pm/[0 III] 88 pm 
line ratio does not show a trend with the ratio [O III]88pm/FIR ratio as it would 
in either of these cases and all of the measured [O III] line ratios fall within the 
range 0.6 - 1.2, consistent with electron densities between 100 - 500 cm -3 . The se- 
quence does show a trend toward lines with lower excitation potentials in the ratios 
[N III]57pm/[N II]122pm (Figure 3) and [O III]52pm/[N III]57pm (Figure 4). No 
FIR fine-structure line emission from species with excitation potentials greater than 
13.6 eV were detected in Arp 220 or in Mkn 231 (Harvey et al. 1999), the two ULIRGs 
for which full ISO LWS spectra were taken. Voit (1992) discussed the possibility that 
even the mid- and far-infrared fine-structure lines would be weak in ULIRGs if they 
are formed in high ionization parameter regions. In such regions with high ratios 
of ionizing photon to electron densities, UV photons are preferentially absorbed by 
dust rather than gas, due to the high column densities of ionized gas in such regions. 
Bottorff et al. (1998) found that the L/^/H/j ratios in such dust-bounded nebulae are 
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Fig. 1. — The full ISO Long Wavelength Spectrometer spectra of six IR-bright galax- 
ies. The spectra have been shifted and ordered vertically according to the equivalent 
width of the [O III] 88 /iin line. The excitation potential, the energy required to create 
the species, is given in eV at the bottom of the figure. From Fischer et al. (1999). 
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Fig. 2. — The [O III]52/im/[0 III]88/im line ratio versus the [O IlljSS/rm line to 
integrated far-infrared continuum flux ratio for the sample galaxies. The dashed and 
dotted lines show the [O III] line ratio in the low density limit (< 100 cm -3 ) and for 
an electron density of 500 cm -3 , respectively (Fischer et al. 1999). 
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Fig. 3. — As in Figure 2 for the [N III] 57 gm/[N II]122/im line ratio. 
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Fig. 4. — As in Figure 2 for the [O III]52 /xm/[N III] 57 gm line ratio. 
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greater than 100 for ionization parameters greater than 10 2 , density = 100 cm 3 , 
and stellar temperatures between 30,000 - 50,000 K. 

Here we present photoionization models with starburst and power law ionization 
sources to predict the strengths of the fine-structure lines in dust-bounded nebulae 
and to compare them with the LWS spectra of the infrared-bright, galaxies. Due to 
the weakness of the infrared fine-structure lines from photoionized gas in ULIRGs, 
ISO mid- and far-infrared spectroscopy produced mostly upper limits. Comparison 
of photoionization models of dust-bounded nebulae with spectra from future space 
missions such as SIRTF and Herschel will help to determine the conditions in the 
photoionized regions of these galaxies. The understanding yielded by these missions 
can then be used to probe the ionized media in high redshift galaxies by the missions 
being discussed in this workshop. 


2. Photoionization modeling 

The photoionization modeling was done using CLOUDY 94.01 (Ferland et al. 
1998) for central power law (Figure 5) and instantaneous starburst (Figure 6) ioniza- 
tion sources. For power law models, the “table” power law option in CLOUDY was 
used. This option produces a continuum with f„ oc v a that is well behaved at both 
high and low energy limits (10 -8 - 10 8 Rydbergs in CLOUDY). An index a = -1.5 
was used for the mid-range (10 /an - 50 keV), while the default indices of +2.5 and 
-2 were used for the low and high ranges, respectively. For the starburst models, we 
used the instantaneous, Salpeter IMF, 3 and 5 Myr aged burst models of Leitherer 
et al. (1999) with solar metallicity and standard mass loss. H II region abundances 
were used in CLOUDY. 

With a central ionizing source, the ionization parameter, defined as the ratio of 
ionizing photons to hydrogen atoms at the inner face of the cloud, U = Qj\. r Kr 1 nc = 
Nh ii chb / c, where Q is the central Lyman continuum rate, n is the density, r is 
the inner radius of the cloud, and as is the Menzel case B recombination coefficient. 
It was varied from 10~' ! — 10 1 by setting Q = 4.5 x 10 54 sec -1 and varying the inner 
radius from 30 - 1600 pc. 

The power law models predict a decrease in the [N III] 57 /rm/[N II] 122 /jin and 
[Ne V]14 fim /[Ne III] 15 fim line ratios in the high U cases. For the starburst models, 
only the far-infrared line ratio decreased at high U, over the range of parameters 
explored. In both cases the line-to-luminosity ratios drop at high U, as expected. As 
Figures 5 & 6 show, high densities also produce both lower values of the line ratios 
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Fig. 5. — Far-infrared (left) and mid-infrared (right) line ratios vs. line-to-luminosity 
ratios for an AGN-like power law ionizing source for moderate (500 cm -3 ) and high 
(10 5 cm -3 ) densities are plotted for ionization parameters U, for log U = -3, -2, -1, 
0, and 0.5 (as shown by the greyscale bar, right). See text for the details of the 
modeling. 
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Fig. 6. — Same as Figure 5 (except for log scale in the line ratios), for ionization by 
instantaneous starbursts of ages 3 and 5 Myrs. 
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and lower line-to-luminosity ratios. Older starbursts can also produce these effects, 
but are unlikely to power ULIRGs (Satyapal et al. 2002). 


3. The effects of high U on photodissociation region diagnostics 

The CLOUDY models presented here show that under high U conditions, line 
emission ratios from ionization states below the dominant stage of ionization at the 
illuminated face of the cloud, i.e. , from constituents deeper in the cloud, are often 
inverted compared with moderate ionization parameters due to the effects of the 
strong grain absorption. Some of the other infrared characteristics of ULIRGs are 
warm 60/100 /im colors, low aromatic feature-to- luminosity and low [C II] 158 /.mi 
line-to-luminosity ratios, although their aromatic feature-to- [C II] ratios are normal 
(Dudley et al. 2002; Luhman et al. 2002). These characteristics may be the result 
of the effects of grains in ionized regions. For example in our own galaxy, Boulanger 
et al. (1988) find that 60/100 /mi colors are higher and small grain emission lower, as 
traced by the IRAS 12 gm flux, within the California nebula H II region than outside 
of it. 


SIRTF and Herschel studies will greatly improve our understanding of these 
diagnostics. SAFIR and future space-borne FIR/submillimeter interferometry will 
be able to exploit them in the high redshift universe to futher our understanding of 
the evolution of galaxies. 


It is a pleasure to thank Gary Ferland for helpful discussions. This work was 
supported by the NASA Long Term Space Astrophysics program and the Office of 
Naval Research. 
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